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, This is a review of strong coupling approaches to grasp the nature of the phase transition 

Qin finite temperature and density QCD. We commence with classics of the center symmetry 
and the Polyakov loop in pure gauge theories. The effective action derived in the strong 
\ coupling limit gives qualitatively plausible behavior for the deconfinement transition at fi- 

nite temperature. We can apply the strong coupling analysis to describe the chiral phase 
7-H ' transition with the help of large dimensional expansion. In order to make this article as self- 

^ ' contained as possible, we elucidate the computational procedure and its physical meaning 

["""»» , in detail. Also the relation between the Polyakov loop behavior and the chiral dynamics, 

. the phase structure of two-color QCD, and future possibilities to be pursued in the strong 

f^) ' coupling analysis are discussed. 

<n : 

£T) . §1. Objectives 

O ; 

■ The aim we are trying to achieve here is to understand the phase structure of 

Quantum Chromodynamics (QCD) in a medium at finite temperature and baryon 
Q_i| density. It is known that QCD would undergo phase transitions between differ- 

ent states of matter, namely, the hadronic phase, the quark-gluon plasma (color- 
deconfined or chiral symmetry restored phase), 1 ) the color-super conductivity or 
.£h ! diquark super-fluidity, 2 ) and some other exotic states like pion and kaon condensa- 

^ ■ tions, ferromagnetic spin alignment, etc. 

The lattice QCD simulation is one of the most prosperous implements to draw 
the right answer among various possibilities. 3 ) The phase structure of QCD or 
QCD-like theories have been intensely investigated by the lattice simulation at fi- 
nite temperature. In many cases the physical quantities of interest in the lattice 
calculation are the Polyakov loop and the chiral condensate. The Polyakov loop is 
supposed to be an order parameter characterizing the deconfinement transition, 4 )' 5 ) 
as explained later. The chiral condensate is a familiar quantity serving as an order 
parameter for the chiral phase transition. 

Even though we could extract the right answer from the lattice QCD results, we 
cannot reach a deep understanding of the QCD phase structure until we correctly 
comprehend, interpret and explain the numerical outputs. In a sense the lattice 
results may be regarded as experimental facts from which we should draw the physical 
meaning. For this purpose the effective model approach would be useful to abstract 
the essence. Also, on a technical level, model studies are indispensable particularly 
at finite baryon density because the lattice simulation cannot avoid suffering from 
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the notorious negative sign problem of the Dirac determinant. 

There are many successful effective models so far, such as the Nambu-Jona- 
Lasinio model, 6 -* the linear sigma model, 7 ) the chiral random matrix model, 8 -* the 
instanton liquid model, 9 ) the Polyakov loop model, 10 ) and so on. In this review, 
we shall draw attention to another effective model description based on the strong 
coupling expansion. Even though the effective action given at strong coupling has no 
guarantee to lead to quantitative agreement with physical observables in the weak 
coupling regime, the strong coupling approach has advantages over other models 
as follows: 1) Because the effective model derived in the strong coupling expansion 
has an obvious connection to the fundamental theory, the result could be exact, 
in principle, in the limit of strong coupling. 2) The strong coupling expansion is 
the most natural approach on the lattice. Since it is formulated on the lattice, 
the correspondence to the lattice QCD simulation is transparent. 3) The strong 
coupling analysis is almost a unique method in which the chiral and the Polyakov 
loop dynamics are taken into account on equal footing. 4) The model parameters are 
few, namely, the lattice spacing a and the current quark mass m q . To our surprise, 
we can reach qualitative and even quantitative agreement with physical observables 
with these two parameters. 

It should be noted here that we must 
be cautious about the continuity between 
the strong and weak coupling regions. 
If there is a phase transition at a cer- 
tain value of the coupling constant, the 
strong coupling physics would not give 
any information on the weak coupling 
physics in reality. Although it is not re- 
solved whether the strong coupling world 
is smoothly connected to the real world 
without any singularity, the string ten- 
sion of the SU(2) gauge theory calculated 
up to twelfth order of the strong coupling 
expansion 11 ^ suggests that physics con- 
tinuously changes from the strong cou- 
pling limit toward the scaling limit in the 
perturbative regime (see Fig. 1). 

§2. Deconfinement and Chiral Restoration 

The deconfinement phase transition is well-defined only in the pure gauge theory 
(gauge theory without quarks in the fundamental representation in color space). In 
this case we can determine whether the vacuum is a confined one or not by putting a 
test static color-charge. If the free energy gain F q diverges, any color-charge isolation 
is not allowed in this system so that confinement is realized. [If a theory contains 
dynamical quarks, F q is always finite due to string breaking effects. It is difficult 
then to characterize confinement by thermodynamic properties. 12 )] 




Fig. 1. String tension in unit of the lattice 
spacing calculated up to twelfth order. 11 - 1 
Two dashed curves represent the extreme 
cases of strong and weak coupling limits. 
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F q can be calculated from the partition function with a test charge carried by a 
static quark q a (f), which can be explicitly written down as 5 )' 13 ) 

F q (r) = -Tin / VA, J^iA; q a (f)\ e'^VlA; q a (f)) (2-1) 

J a=l 

in the Weyl gauge (^4o = 0). V is the projection operator onto the physical state 
which should satisfy the Gauss law, that is; 

-ign q ^- q yA f ,qP(r)) = +g^-6^6(x - r)||A) = (2-2) 

with the canonical momenta Fl % £ and Fl q for the gauge and quark fields respectively. 
Di stands for the covariant derivative. The Gauss law constraint needs a Lagrange 
multiplier which can be eventually regarded as A4 (temporal component of the gauge 
field in Euclidean space-time). 14 ) Then in the imaginary-time (Matsubara) formalism 
the free energy can be expressed as 

F q (r) = -TlnJvA^e- s[A] Tr c Texp ig dx 4 A 4 (r, 24) = -Tln(Tr c L(r)), (2-3) 

where T is time ordering and Tr c is the trace in color space. The Polyakov loop is de- 
noted by L. If we introduce a pair of test charge and anti-charge in the same way, we 
can attain the inter-quark potential, F qq , expressed by the Polyakov loop correlation 
function. In Table I the basics on the Polyakov loop behavior are summarized. 





confined (disordered) phase 


deconfined (ordered) phase 


free energy 


F q — 00 


F q < 00 






F<aW ~ -e- m ^/[r\ 


Polyakov loop 


(Tr c L) = 


(Tr c L) / 




(Tr c tf(0)Tr c L(r)) ^0 


(Tr c L t (0)Tr c L(r)) (|f| -> 00) 



Table 1. Confined and deconfined phases in the pure gauge theory. 

As is clear in Table I, the expectation value of the Polyakov loop serves as an 
order parameter to identify the deconfinement transition. In many cases a phase 
transition is linked with the spontaneous symmetry breaking. As a matter of fact, 
the Polyakov loop behavior is prescribed by the center symmetry , 5 )> 13 )> 14 ) The cen- 
ter symmetry in the SU(iV c ) gauge theory is defined by the gauge transformation 
satisfying a boundary condition twisted by an element, z, of the center Z(A r c ). 15 ) 
The point is that such a gauge transformation does not affect the periodic boundary 
condition for the gauge field, while the Polyakov loop picks up the element z. As a 
result the expectation value of the Polyakov loop vanishes if the vacuum preserves 
the center symmetry. Hence, concerning the deconfinement transition, the SU(iV c ) 
gauge theory can be regarded as a spin system with the global Z(A r c ) symmetry 
described in terms of the Polyakov loop "spin" variables. The universality argu- 
ment tells us the critical properties of the deconfinement transition which have been 
confirmed in the lattice simulations. 3 ^' 13 )' 16 )' 17 ) 
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One might have thought that it should be curious that the disordered (ordered) 
phase lies in the lower (higher) temperature side. 18 ) This seeming contradiction 
comes from the duality transformation connecting the gauge theory to a spin system. 
Going back to the original paper written by Polyakov, 4 ) we can find the explicit 
transformation for the partition function (notation is slightly changed); 



T 

2g 

(2-4) 

The left-hand-side is the partition function given by Hamiltonian in the strong cou- 
pling limit on the lattice. The electric flux string, n I(i , dominates over the magnetic 
fluctuation at strong coupling. The Gauss law is implemented by the Lagrange 
multiplier <p (that is A4 in Eq. (2-3)). The right-hand-side is derived by Poisson's 
resummation formula (duality transformation) up to an overall factor and is nothing 
but an approximated form of the XY spin model. From the above expression it can 
be clearly understood why the effective "spin" theory of the Polyakov loop has the 
inverse temperature. 

In the rest of this section we shall briefly summarize the definitions inherent in 
the lattice description. The Wilson action of the gauge field is given by 19 ) 

S gP] = E I 1 " ^ReTr c CV(*)} (2-5) 

with the plaquettes; U liU (x) = ut(x)Uj l (x+u)U u (x+fi)U ll (x). The coupling constant, 
g 2 , only appears as an overall coefficient, from which the strong coupling expansion 
is derived. The Polyakov loop can be equivalently written on the lattice as 

N T a 

L(x)= H U d (x,x d ) (2-6) 

x d =a 

in d— dimensional space-time. N T is the temporal extent and related to the temper- 
ature by T = l/N T a with the lattice spacing a. 

As for the chiral symmetry, the lattice fermion suffers from the notorious species 
doubling problem. Among ideas to avoid having 2 rf -fold doublers, we shall limit our 
discussion in this review to the staggered formalism (Kogut-Susskind fermion 20 )) that 
has advantage in providing a simple description of chiral symmetry. The staggered 
fermion x( x ) is derived from the Dirac fermion ip{x) by the following transformation; 

1>(x) = r(x) x (x), i>{x) = x{x)r\ x ) : r(x) = 7 r72 2 ---7?- (2-7) 

Then Dirac's gamma matrices are reduced to be unity and we can forget about the 
Dirac indices. As a result 2 d -fold doublers are diminished by 2^ d / 2 \ The remains 
are usually interpreted as flavor degrees of freedom. Since the flavor number might 
be sometimes confusing as it is, we will write nf to mean the original flavor num- 
ber associated with staggered fermions and Nf to mean the flavor number in the 
continuum limit (i.e., Nf = 4nf for d = 4). 



Strong coupling analysis 



5 



It should be noted that the chiral symmetry is given by a U(nf)xU(rtf) rotation 
on alternate lattice sites. Although the symmetry pattern looks different, the chiral 
condensate, (xx): serves as an order parameter for the chiral symmetry breaking; 
U(nf) x U(nf) — > Uv(wf). The action for the staggered fermion is given by 

Sf[U,x,x] =m q ^2x(x)x(x) + ^ ^ 7] ll (x)x(x){U ll {x)x(x+fi)-Ul l (x-fi)x(x-fi)} 

x x,[i>0 

(2-8) 

with rju(x) = (— i s ) x i+ x 2+— x n-i_ The price for reducing spin indices is that the flavor 
contents are scattered and mixed on the adjacent lattice sites. 



§3. Deconfinement Transition 



The lattice gauge theory was originally formulated by Wilson 19 ) with the aim 
of giving a plain explanation of color-confinement in the strong coupling limit. The 
strong coupling expansion on the lattice is achieved by statistical analog at high 
temperature. The formulation has been developed systematically. In this article we 
will not go into mathematical and technical details in general. If necessary, one can 
consult an exhaustive review written by Drouffe and Zuber. 21 ) 

In this section we shall focus on the strong 
coupling study on the deconfinement transition 
for which the order parameter is given by the 
Polyakov loop. We can attain the effective action 
in terms of the Polyakov loop by integrating over 
other degrees of freedom. As seen from (2-5), the 
large g 2 expansion is systematically available by 
the Taylor expansion of exponential of Sq . In the 
leading order of the strong coupling expansion 
the partition function with the Polyakov loop (or 
Ud) left unintegrated is calculated as (see Fig. 2) 




Fig. 2. Leading-order contribution 
to the effective action with re- 
spect to the Polyakov loop. The 
cylinder represents temporal S 



and spatial R manifold. 



N T a 



1+ E ^ II U d (x,x d )Tr c H U d (x+ja,x' d ) 

x,j<d x d= a x' d =a 



1 \N. 



9 2 N C 



+ (h.c.) + • 
(3d) 



Note that 1/N C emerges as a result of the group integration. This expression imme- 
diately leads to the effective action as 



S$[L] = -lnZ^ } = -e"™/ T ^Tr c L(x)Tr c Lt(y), 



r(G) 



(3-2) 



where n.n. is the abbreviation to mean the nearest-neighbor interaction. We have 
used T = l/N T a and the familiar expression of the string tension in the strong 
coupling limit; era 2 = ln[g 2 N c ]. 
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The effective action (3-2) describes the (d— l)-dimensional classical spin system 
with the spin variable represented by Tr c L and the exchange interaction by e - aa / T , 22 ) 
Although (3-2) looks quite simple, non-trivial complexity can originate from the 
matrix nature of the Polyakov loop. In the simplest mean- field analysis the matrix 
nature is taken into account by the Haar measure explicitly, that is the Jacobian 
associated with the variable transformation from At to L. 23 ) Then the effective 
potential can be approximated as 

V^\L]/T = -2(d - l)e-™/ T |Tr c L| 2 - hi Mm*r[L], (3-3) 

where 

, tJ .. flnfl - |Tr c L/2| 2 l for iV c = 2, , N 

lnM H aar^=< r o « a, (3-4) 

1 J \ln[l-6|Tr c L/3| 2 +8Re(Tr c L/3) 3 -3|Tr c L/3| 4 ] for iV c = 3. V ' 

The effective potential for N c = 2 gives a second-order phase transition, 24 ) while the 
Z(3) center symmetry allows a cubic term so that the phase transition is of first- 
order for iV c = 3. 25 ) Interestingly enough, it is obvious from (3-3) that confinement 
as well as the Z(iV c ) nature stem only from the Haar measure, as is consistent with 
the arguments that the Haar measure is essentially important for non-perturbative 
aspects of QCD. 26 ) 

Unfortunately straightforward generalization to iV c > 3 cases is not available. 
Instead, as usually performed in the spin system, the Weiss approximation works 
well for arbitrary N c . In this approximation, we take account of the fluctuation of 
the individual Polyakov loop surrounded by a constant mean field. 21 )' 22 )> 27 )- 29 ) This 
results in a first-order phase transition for N c > 3. In particular we can find the 
critical coupling 2(d— l)e -CTa / Td = 0.806 for iV c = 3. Once empirical values, the string 
tension a = (425 MeV) 2 and the lattice spacing a -1 = 440 MeV (see next section), 
are substituted, the deconfinement temperature is estimated as = 204 MeV. One 
can, otherwise, fit the empirical value, = 270 MeV, by choosing a -1 = 333 MeV. 29 ) 

In short summary, the strong coupling expansion yields a simple form of the 
effective potential in terms of the Polyakov loop. It leads to even quantitatively 
acceptable results as well as qualitative agreement with gross features of the decon- 
finement transition. Since it is easy to handle, the effective potential (3-3) or more 
sophisticated one 28 )' 29 ) can also be a useful ingredient to build an effective model 
and probe underlying physics in model studies. 30 ) 



§4. Chiral Phase Transition 



4.1. Spontaneous chiral symmetry breaking at T = /i = 

Now we shall deal with the system with dynamical quarks. The center symmetry 
is lost and the deconfinement transition is obscured then. Instead, the chiral sym- 
metry plays an important role. In this subsection we will describe the spontaneous 
chiral symmetry breaking at T = fi = mainly according to the formalism given by 
Kawamoto and Smit 31 ) and Kluberg-Stern, Morel and Petersson. 32 ) 
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$ ♦ I 
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M(x) M(x+\i) Bfy) B(x+\i) 

Fig. 3. Intuitive graphical representation of the C7 M integration. 

In the g 2 —* oo limit the gauge action can be neglected. We can accomplish 
the integration with respect to on each lattice site by expanding the exponential 
of (2-8). Such series of the Taylor expansion corresponds to the large dimensional 
expansion, as explained below. For our present purpose to explain how it works, 
we only keep the leading order contributions. Then the XJn integration picks up the 
mesonic and baryonic (color-singlet) composites as (see Fig. 3) 



dUfj,(x) exp 
1 + 2^ Tr f M {x)M{x + £) + ••■ 

f l-\No(N c -l)/2^ 1l( x ) ^ N ' 



Tr { {B(x+p,)B(x) + (-l) N °B(x)B(x+ji)} + • • • , (4-1) 



where the mesonic and baryonic composites with flavor indices are defined by 



MM 



1 id Nc 

j\T V 9 Xi X 3 



B 



a=l 
d N «/ 4 



(4-2) 



£ai_a2---a Nc Xi 1 Xi 2 * * * A 



(x). 



oi,...,«jv 



Here we can understand where the 1/d expansion emerges from. The large 
d limit would make sense under the condition that the mesonic propagation has 
a finite amplitude. Since the index fi runs over d directions, the hopping term 
M(x)M(x + ft) should be divided by d. In other words, the quark fields, x an d X, 
should be normalized by d 1 / 4 like (4-2). Consequently, the Taylor expansion in terms 
of x an d X leads to the systematic 1/d expansion. 

Although the baryonic term may be important when the system has a non- 
vanishing baryon chemical potential, we shall drop it from (4-1) and examine the 
spontaneous breaking of chiral symmetry. [Note that the baryonic term is of higher 
order in the 1/d expansion for N c > 3.] In order to transform variables from the 
quark fields x to the meson field A, we linearize the four-fermionic interaction by 
inserting an auxiliary field in the following way; 



exp 



N c TrA~M(x)V M (x, y)M(y) + 2m q M{x)5 3 

N c J2TrAh(x)V M (x,y)X(y) - (X(x)V M (x,y) + 2m q 6 X)y )M{y) 



VX exp 



(4-3) 
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where Vm(x,u) = (l/2cZ) Y^^x,y+fi + $x,y-[i) represents the mesonic hopping prop- 
agation between the nearest neighbor sites and m q = m q / y/2d. 

Now that the effective action is linearized in regard to the composite field M(x), 
the integration with respect to the quark fields, x( x ), is straightforward to be carried 
out resulting in the effective action in terms of A(x); 

Srf W = ^ E Tr f X ^ X ) V M (*> V)KV) "iVcE Tr f ln [E Kv)V{y, x) + 2fh q ] . (4-4) 

x,y x y 

Here A can be regarded as the meson field and its condensate, which is denoted by 
<5 y A, generates an additional contribution to the quark mass. The stationary value 
of A is fixed by the extremal condition of the effective potential, 

Kff } W = ^pA 2 - N c n { In [A + 2m q ] , (4-5) 

to lead to a finite scalar condensate, A = — fh q + yjl + to 2 ,, which brings about the 

spontaneous breaking of chiral symmetry. Actually A is directly related to the chiral 
condensate, 

<*<*>*<*» = -]^-dW PAe ~* w ' (4 ' 6> 

by a simple proportionality relation; \P q = —X^/2/d (see a reference 29 -* for details). 
We should remark that the logarithmic term in the effective potential (4-5) makes 
infinite barrier around the symmetric vacuum (A = rh q = 0) and this logarithmic 
singularity alone is responsible for the symmetry breaking in the present approach. 

In the next place let us look into the meson spectrum acquired in this framework. 
The kinetic term of the meson fluctuation, 5X(x), around the stationary condensate 
provides us with the propagator of meson excitations, whose poles correspond to the 
dispersion relations physical particles should satisfy. From the effective action (4-4) 
we have the kinetic term as 

-^.^Tr t 8X(x){v M (x,y)+ - + ^_ 2 V^x,y)}5X(y). (4-7) 

x,y Q 

To derive the information on the meson spectrum we transform the meson propagator 
into the representation in the momentum space. Then the inverse propagator is 

a ^2fh) 2 £ eikXVM ^ y^~ X + *™ q )%0 + V M (y, 0)} 

^ 9' x,y 

oc V M (k){(X + 2m q ) 2 + V M {k)} oc - E cos k v + (A + 2m,) 2 , (4-8) 

where we used the Fourier transformed hopping propagator in the momentum space, 
VAf(fe), which is obtained as 

V M (k) = E e ikx V M (x, 0) = - d E cos k„. (4-9) 

X fj, 
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Now that the meson propagator in the momentum space is known, we can read the 
meson masses putting the momentum k^ = (0, iM) + tt5^ into the propagator. The 
last term <5„ takes either 1 or representing the fermion doublers which are absorbed 
into the flavor degrees of freedom. From the condition that the inverse of the meson 
propagator has zeros, we have the meson spectrum in the present approach as 

cosh M p = d{{\ + 2rh q f - 1} + 2p + 1. (p = 0, 1, . . . , d - 1) (4-10) 

Mo and Mi are identified as the masses of the lightest state belonging to 0~ (pion) 
and mainly 1~ (p meson), respectively. M2 and M3 are regarded as 1 + (cii meson) 
and + (ao meson) with significant mixtures with 1~ and 0~. 33 ) In order to reproduce 
the physical values, = 140 MeV and m p = 780 MeV *\ the model parameters 
are fixed as 

m, = 8MeV, aT 1 = 440 MeV, (4-11) 

which give other physical quantities as listed in Table II. We would like to conclude 
this subsection with quoting from Kluberg-Stern et al. - "We believe that this is not 
a bad result owing to the crudeness of the model. " 





9*=0 


Physical values 


M 2 


910 MeV 


m ai = 1260 MeV 


M 3 


1010 MeV 


m ao = 980 MeV 


M B 


1380 MeV 


m N = 940 MeV 


% 


-(390 MeV) 3 


-(240 MeV) 3 


u 


190 MeV 


93 MeV 



Table II. Hadron spectrum inferred from 
the strong coupling analysis. The 
baryon mass, Mb, can be obtained by 
inserting the baryonic auxiliary field to 
(4-1) and f n from the PCAC relation. 



4.2. Chiral restoration at T / and n / 

Since the 1/d expansion works quite well to describe the spontaneous chiral sym- 
metry breaking at zero temperature, it is interesting to see what kind of phase struc- 
ture appears at finite temperature and baryon density. The extension to the evalu- 
ation in a hot and dense medium is straightforward as found in literatures. 28 )' 29 )> 34 ) 
The temperature dependence originates from the boundary condition in the tem- 
poral direction and the baryon chemical potential dependence is introduced by the 
replacement; Ud — > e^Ud (// being the quark chemical potential). 35 -* These facts 
suggest that we should carefully treat the integration along the temporal direction 
apart form the others in the spatial (d — l)-dimensional coordinates. We have two 
strategies to embody it. 

One is a physically preferable but technically difficult method. First, all Un 
(including Ud) are integrated by means of the 1/d expansion. The resultant action 
is (4-1) with the baryonic term in the temporal direction multiplied by e^^. Then, 
after inserting the mesonic and baryonic auxiliary fields (A and b), the integrations 
with respect to x an d b are performed to obtain the effective potential in terms of 
A. This computational procedure works well at T = and /i / 0. The physical 
intuition that the \i dependence comes from the baryonic excitations would make 

We adopt not m p — 770 MeV but this value faithfully according to Kluberg-Stern et al. 32 ^ 
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sense. Unfortunately for the T / system, however, we have to take account of not 
only the baryonic but also the mesonic loops explicitly and the calculation becomes 
too complicated. 

The other is much simpler to handle. First, Ud are left unintegrated and the 
1/d expansion is applied only for the U{ integration. Then, after inserting A, the 
integration with respect to x anci Ud are exactly carried out. It should be mentioned 
that in this method the baryonic term in (4-1) is not necessary to deal with the \i 
dependence. In principle, the x integration with unintegrated Ud contains in itself 
the thermal excitation of both mesonic and baryonic composite states. The problem 
is that the physical meaning is not transparent especially in the confined phase 
because the thermal excitation is given by quarks in this method. 

Here we will explain the latter method in detail neglecting the baryonic hopping 
term. The integration over Ui(x) is done by the same procedures as in the case at 
zero temperature. Then the partition function can be written as 



VU d VxVx exp 



l -Y,^)x{x){^U d {x) X {x + d) - e-»u\{x-d) X {x-d)} 

(4-12) 



2 

X 



^ X(x)V M (x, y)\{y) + N c ^ (\(x)V M (x, y) + 2m q 8 x , y ) M{y) 



2 

x,y x,y 

where VM(x,y) represents the hopping propagator only in the spatial directions. In 
the mean-field approximation the meson field, A(x), is simply replaced by constant 
S^X. Then the functional integration with respect to the quark field, x( x ), is reduced 
into the one-dimensional problem and it can be easily manipulated. After the x 
integration, the expression is simplified as 

r Nc 

/ VU d ]J [ 2 cosh(N T E) + 2 cos{9 a - i/x) "* (4-13) 

a=l 

in the Polyakov loop gauge (U d = diag(e l6)l , . . . ,e 10N ^)). For simplicity we will focus 
on the rif = 1 case (iVf = 4) from now on. Then the integration with respect to x 
results in the effective potential; 



, /AT AT \ sinh[(iV c + 1)N T E] 
2 cosh iV T iV c/ u + U * AT ' T J 



(4-14) 



with the definition; E = smhT 1 ^^ (d — l)/2 + m q ). 

Now we can draw the phase diagram derived from (4-14) with m q = (chiral 
limit) and iV c = 3. As long as \x is not large, the chiral phase transition is of second- 
order. The phase boundary can be determined by the condition that the potential 
curvature around A = (denoted by Co) passes across zero. When /j, gets larger, 
the phase transition becomes of first-order. As a result we reach the phase structure 
depicted in Fig. 4. What seems interesting is that in some region of fx the chiral 
phase transition undergoes twice with increasing temperature unlike the results in 
other model studies. 
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c„<o 



Ist-order 



0.2 0.4 0.6 

Quark Chemical Potential // 



Fig. 4. Phase diagram derived from (4-14) 



with m a 



and N c 



3 in unit of 



the lattice spacing. Co is the potential 
curvature at A = 0. Thick dashed and 
solid curves represent the second and 
first order boundaries respectively. 



§5. Relation between Two Phase Transitions 



In the last section we have seen that the Ud integration is dealt with separately. 
Since the Polyakov loop consists of Ud, it is easy to put the Polyakov loop into the 
formalism. Once we construct an effective model with the chiral order parameter 
and the Polyakov loop, we can investigate the relation between the deconfinement 
and chiral phase transitions. An interesting question is why these two transitions 
have been observed at the same temperature in the lattice QCD simulations. 17 )' 36 ) 

The effect of dynamical quarks on the deconfinement phase transition was first 
studied by Green and Karsch 37 ) with the hopping parameter (heavy quark) expan- 
sion. Chiral symmetry was seriously taken into account by Gocksch and Ogilvie, 28 ) 
followed by the finite density extension by Ilgenfritz and Kripfganz. 38 ) TheGocksch- 
Ogilvie model is defined by the effective action; 29 ) 



5 eff [L,A] = -e-™/ T ^Tr c L(x)Tr c Lt(y) 

n.n. 

+ H*)Vm(x, y)X(y) - ^ £ Tr c In 

x,y x 



cosh(N T E) + -(L + ]J) 



(5-1) 



This can be understood as a combination of (3-2) and (4-14). The effective potential 
V^ff [Z, A] (I being the expectation value of the Polyakov loop) obtained from (5-1) 
in the mean-field approximation has an interesting property, that is, V c s[l = 0, A] 
inevitably has instability at A = and thus it leads to a non-vanishing chiral con- 
densate. Roughly speaking, the confined vacuum breaks chiral symmetry at any tem- 
perature. As emphasized in a chiral effective model study with the Polyakov loop, 30 ) 
this property suggests T x > Td (the chiral restoration temperature is higher than the 
deconfinement temperature). Remembering that Td ~ 270 MeV without dynamical 
quarks is higher than T x ~ 150 MeV, we can expect to have Td = T x in the presence 
of dynamical quarks. As a matter of fact, we can see that in the Gocksch-Ogilvie 
model the simultaneous crossovers are attributed to this mechanism, as shown in the 
left of Fig. 5. 
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Fig. 5. The Polyakov loop and the chiral order parameter as functions of the temperature calculated 
in the Gocksch-Ogilvie model. 29 ' The left figure is for the system with dynamical quarks in the 
fundamental representation in color space. The right figure is for the system with adjoint quarks. 
Thick curves are for the physical value of the quark mass and thin solid and dashed curves are 
for m q = 1 MeV and m q — 50 MeV respectively. 



§6. QCD-like Theories 

6.1. QCD with Adjoint Quarks (aQCD) 

Adjoint quarks have the same color structure as gluons and thus never break the 
center symmetry. In the lattice simulation of QCD with adjoint quarks (aQCD), as 
expected, the deconfinement transition is found to be of first-order. 39 -* The Gocksch- 
Ogilvie model has been generalized to include adjoint quarks and the model study 
gives a satisfactory explanation on all the qualitative features. 29 ^ The right of Fig. 5 
is the result in the Gocksch-Ogilvie model with adjoint quarks. 

6.2. Two-Color QCD 

Two-color QCD (QCD with iV c = 2) has special properties. As explained before, 
the staggered fermion has the global U(nf) xU(nf) symmetry for m q = 0. For iV c = 2, 
the chiral symmetry is graded to U(2nf). 32 )' 40 ) This is because the color SU(2) group 
is pseudo-real and the action possesses Pauli-Giirsey's symmetry, 41 ) which also makes 
possible the lattice simulation at finite density. 

Introduction of a finite chemical potential /j, or quark mass m q explicitly breaks 
the U(2rif) symmetry. Listed in Table III are the symmetries realized in various 
circumstances and their breaking patterns for nf = 1. 





m = 


m j= 


n = 


U(2) 

broken to U(l) with 3 NG modes 


U V (1) 
not broken 




U v (l) x U A (1) 
totally broken with 2 NG modes 


U V (1) 

totally broken with 1 NG mode 



Table III. Symmetry realized in the single-component staggered-fermion action for iV c = 2. Possible 
symmetry breaking patterns and the number of Nambu-Goldstone (NG) modes are also listed. 
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To our surprise, the chiral condensate, (xx)> °f two-color QCD vanishes in the 
chiral limit {m q = 0). Instead, the diquark condensate, {xx)i pl avs an important role. 
Note that XX = e ai3X°X 13 is a gauge invariant operator for N c = 2 and it corresponds 
to the baryonic contribution in (4-1). The phase structure was first investigated in 
the strong coupling approach by Dagotto, Karsch and Moreo. 42 ) The analytical and 
numerical studies have been inexhaustively explored within the framework of the 
strong coupling limit. 43 )' 44 ) The results are consistent with the lattice simulation. 45 ) 
The order parameters and the phase diagram are shown in Fig. 6. 




M 



Fig. 6. Left: Chiral and diquark condensates (a and A respectively) for m q = 0.02. Dotted curve 
represents the magnitude of condensation; \/a 2 + A 2 . The mixing angle between a and A is 
controlled by m q . If m q — then a = and the condensation is saturated by A. Right: 
Phase boundaries in the 3-dimensional (m q ,T,[i) space. Below the critical surface, A takes a 
non-vanishing value. 44 - ) 

§7. Beyond Understanding 

We have looked over various aspects of the strong coupling approach to QCD. 
Since the celebrated paper by Wilson, there are quite a few contributions to clarify 
QCD physics in the strong coupling limit. Nevertheless, there still remain several 
interesting questions to be answered in the strong coupling analysis. We shall enu- 
merate some possibilities here; 

1) Glueballs at finite temperature - In the pure gauge theory, the Polyakov loop 
is an order parameter for the deconfinement transition. However, the physical inter- 
pretation and the dynamical description (real-time evolution) of the Polyakov loop 
are not clear in Minkowski space-time. In principle, we can define the deconfinement 
transition in terms of physical excitations, namely, glueballs. 46 ) At zero temperature, 
glueballs are well described in the strong coupling expansion. 47 ) A finite-temperature 
extension would give us concrete information on electric and magnetic glueballs near 
the deconfinement transition temperature. 48 ) 

2) Relation between the deconfinement and chiral phase transitions at finite 
density - Since the lattice QCD simulation has given no definite answer at finite 
density, the model study preceding future lattice simulations at finite density would 
be especially useful. Because, unlike the temperature, the baryon density effect 
explicitly breaks the center symmetry, it could be expected that the Polyakov loop 
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behavior should be significantly affected at finite density. 

3) Color-super conductivity on the lattice 50 ) - Although the lattice simulation 
is not applicable at finite density, we believe that the sign problem could be gotten 
over someday. Then an interesting and important problem would be how to see the 
color-super conductivity on the lattice. Because of Elitzur's theorem, 49 ) we have to 
make a gauge invariant order parameter or choose a certain gauge fixing condition. 

4) Phase structure of the system at finite isospin density - It is interesting 
to apply the strong coupling approach to investigate the phase transition at finite 
isospin density, 51 -* for there are numerical outputs from the lattice simulation. 52 ) In 
order to incorporate the flavor structure the Wilson fermion is more suitable then. 
The strong coupling analysis with Wilson fermions should be developed. 53 ) 

Beyond understanding the existing lattice data, we would like to emphasize 
that the strong coupling analysis could be a powerful tool to give some guideline or 
■prediction to future progress in the lattice approaches. I hope that this contribution 
presented here will provide some clues to open the way for finite density QCD. 
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